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ABSTRACT 

We propose a novel method for observing the gravitational wave signature of super- 
massive black hole (SMBH) mergers. This method is based on detection of a specific 
type of gravitational waves, namely gravitational wave burst with memory (BWM), 
using pulsar timing. We study the unique signature produced by BWM in anomalous 
pulsar timing residuals. We show that the present day pulsar timing precision allows 
one to detect BWM due to SMBH mergers from distances up to 1 Gpc (for case of equal 
mass 10*^ Mq SMBH). Improvements in precision of pulsar timing together with the 
increase in number of observed pulsars should eventually lead to detection of a BWM 
signal due to SMBH merger, thereby making the proposed technique complementary 
to the capabilities of the planned LISA mission. 

Key words: gravitational waves - galaxies: evolution - (stars:) pulsars: general - 
cosmology: miscellaneous 



INTRODUCTION 



The prospects of detecting Rravitational waves (GWs) in the coming decade are looking ever more promising (iGrishchuk et al] 
l200ll ? ICutIer fc Thornell2002l : ISathvaprakash fc Schutj|2009h . There is currently a considerable experimental effort to detect 
gravitational waves in a wide range of frequencies. At high frequ encies v ~ 10^ — 10'^ Hz the g round based laser interferometric 
detectors, such as the currently operational LIGO (|LIGO||2009| ) and VIRGO (|VIRGOll2009l ). especially in their "Advanced" 
config uration, should be able to observe GWs from a wide range of astrophysical sources. The planned space interferometer 
LISA (|Baker et al.ll2007l ) would be able to observe GWs in the region of frequencies i/ ~ 10 * — 10 ^ Hz. At intermediate 
range of frequencies v ~ 10~* — 10~® pulsar timing measu rements are a ver y strong tool for observing th e GW signature. 
Several pulsar timing array (PTA) projects, such as PPTA (|Manchesteill2007l L EPTA jjanssen et al.ll2008l ') and NANOGrav 
l|Demorest et al.ll200 9l) are collecting data and yielding an unprecedented and every increasing sensi tivity in i ntermediate 
frequency range. Furthermore, the implementation of the planned Square Kilometer Array (SKA) ( Cordes et al.|[20 04;^ radio 
telescope would further improve the sensitivity of pulsar timing measurements to gravitational waves. Finally, at the lowest 
frequencies v ~ 10~^* — 10~^® Hz considerable efforts are being made to measure the impri nts of relic gravitatio nal waves in 
the temperature and polarization anisotropics of the cosmic microwave background (CMB) ([Keating et al.|[2006l ). 

In what follows we s hall consider pulsar timi ng, in particular timing of millisecond pulsar (MSPs), as a tool to study 
gravitational waves (see (|Hobbsl [20051 : [jenet et al.|[2005i ) for a recent overview). A GW passing between the pulsar and the 
observer on Earth will leave a modulating signature on the observed period of the pulsar. For thi s reason, the analysis of 
pulsar timing residuals gives a uniqu e opportunity to observe GWs l|Estabrook fc Wahlguistl [19751 : ISazhinlll978l : iDetweileij 
I1979I : iBertotti. Carr fc Reedll98^ : iKopcikin 1997). A particu larly attractive method for detection is to cros s-correlate the 
timing signal from di fferent pulsars (HcUings fc Downs 1983'). This method led to the proposal for a PTA (|Romanilll989l : 
[Foster fc Backeij|l990l). which is aim e d at o bserving a number of pulsars distributed on the sky at regular time intervals over 
a long time span fsee lVerbiest et al.l (|2009l ) for recent review). 

Current pulsar timing observatio ns are yet to detec t gravitational waves. However, current upper bounds on the stochastic 
< 2 • 10~* ijjenet et al.ll200^ ) (in terms of the ratio of energy density per unit logarithmic frequency 
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interval to the cosmological critical density at frequency f = l/8yrs) already place interesting limits on cosmic string models. 
Moreover, in the near future, it is Hkely that the GW background of coalescing extragalactic sup ermassive black hole binaries 
and possibly the background of relic gravitational waves will be discovered ( Manc hester! 120071 ) . Apart from the stochastic 
backgrounds, pulsar timing also has the potential to observe GWs from in dividual sources. I n particular, there has been 
considerable interest in GWs emitted by massive black hole binaries systems (|jenet et al.ll2003 : ISesana et al.ll2009l ). 

In the present work we shall be interested with indi vidual sources which emit gravitational waves of a particular class , 
namely gravita tional wave bursts with mem ory f BWM) JZel'Dovich fc PolnarevI Il974l : ISmarJ Il977l : iKovacs fc Thornd Il978l : 
iBontz fc Pric3 ll979'; Bragi nskv fc Grishchuklll979l : iBraginskv fc Thornelll987l ). and analyze the prospects for their detection 
in pulsar timing measurements. In general, BWM are characterized by a rise in the gravitational wave field hj^^ from zero, 
followed by oscillatory behaviour for a few cycles and settling dow n at a final non-zero value Ahf^^ after a characteristic time 
6t known as the duration of the burst jBraginskv fc Thorne|[l987l ). The permanent change Ahjj'" is usually referred to as the 
burst's "memory" . BWMs are produced in situations where there is a net change in the time derivatives of multipole moments 
characterizing the system, for example during flyby of massive bodies on hyperbolic trajectories or an asymmetric supernova 
explosion. A BWM where the net change of the moments is attributed to the gravitons released duri ng the burst is known as 
the Christodoulou effect (|Pavndll983l : IChristodouioulll99ll : [Thorndll992l : lBlanchet fc DamouJ[l993 ). Whatever the physical 
reason, the existence of the permanent offset Ahfj'" in BWMs makes them particularly interesting for pulsar timing studies. 
As we shall show below, the burst's memory leads to a characteristic signature in pulsar timing residuals which accumulates 
linearly with time. 

The structure of the paper is as follows. We begin in Section [2] with a simple order-of-magnitude estimation of the 
potential sources of BWM, and compare them with the expected sensitivity of pulsar timing measurements. In this section, 
we study the typical amplitudes and possible event rates for potential signals. The simplistic analysis shows that, in the 
context of pulsar timing, the most observationally promising sources of BWM are extragalactic coalescing supermassive black 
holes binaries. In Section [3l we study the observational prospects for the detection of the BWM signal in significant detail. 
Using a simple analytical model for a BWM signal we calculate the pulsar timing residual. Following this, we evaluate the 
expected signal-to-noise ratio, and the expected number of BWM events that could be observed in a typical pulsar timing 
experiment. Finally, we conclude in Section [4] with an overview and discussion of the main results of this work. 



2 POTENTIAL SOURCES OF GRAVITATIONAL WAVE BURSTS WITH MEMORY 
2.1 Typical strength of BWM signal 

In an astrophysical context, BWM typically occur in burst events which are accompanied by considerable amount of mass 
or radiation ejected in an asymmetric fashion. A simple formula for estimating the characteristic amplitude of BWM 
l|Braginskv fc Thornd[l987l ) reads: 

r V c / 

where rg and v are the Schwarzschild radius and velocity of aspherically ejected parts respectively, and r is the distance to 
the source of the burst. The effect is maximal when the mass is ejected at maximal speed, i.e c, corresponding to ejection of 
photons or gravitons. In the evaluations below we shall work in units [G = c = 1]. 

An obvious source of GW BWM is a core-collapse supernova. In a typical scenario , the value of asymmetri cally radiated 
energy is AiJ^ad ^ W~^Mq (inferred from the velocities of neutron stars, see e.g. iNazin fc Postnovl (| 19971 )). leading to 
an estimate /iiMpc ^ lO"'^^, for the typical BWM strai n fo r a co re- collapse supernova at a distance of 1 Mpc from the 
observer. In fact, these estimates are overoptimistic fsee lOttI (|2009l )). Taking into account the expected event rate of core- 
collapse supernovae of a few per century in a Milky Way type galaxy, one can conclude that the sensitivity of pulsar timing 
measurements will not allow the detection of such a signal (see Section r3.3 | for a discu ssion of the signal-to- noise ratio). 

Another interesting class of burst events are black hole (BH) mergers T Favatall2009l ). Depending on the angular momentum 
of the merging BHs, up to several percent of the total mass of the system can be radiated non- spherically during the burst 
l|Reisswig et al.ll2009l ). The masses of BHs range from stellar to several billion solar masses. However, from ([l]), we conclude 
that the most interesting candidates in the present context are the mergers of super massive black hole (SMBH) binaries, 
with typical masses Msmbh 5S 10* M0. Assuming that 10% of the mass is radiated during the burst, for the BWM strain at a 
distance of 1 Gpc due to a SMBH merger with total mass 10* Mq, we arrive at an estimate /iicpc — 10~^^. As we shall show 

below, this strain is within the reach of pulsar timing measurements. 

This estimate can be made more precise. The BWM amplitude from SMBH merger is given by Eq. (5) from iFavatal 

l|2009l) : 

h = — — sm 61(17 + cos 9), where h = . , • (2) 



3847rr ^ " \ M 

In the above expression, M = Mi + M2 is the total mass of SMBH, 77 — ^^^^^^'''^ and is the angle the binary angular 
momentum and the line of sight. This formula can be rewritten in the form: 

I mem Ai?rad . 2 ^7 , 2 /i\ /o\ 

h = sm 6'(17 + cos 6' . (3) 

24r- 
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Averaging over 9 yields: 

^ ,mom ^ 69 Ai?rad ^ A_Brad 

Estimates of AE^^d can be obtained from calculations of lReisswig et all (|2009h which give Ai?rad from 3.6% to 10% M. In the 
present work, we shall use the mean value of AiSrad ~ 7-10"'^ M. For equal-mass SMBH system with Mi = M2 = m = 10* Mq) 
and zero orbital eccentricity at a distance of r = 1 Gpc, we obtain: 

'■-=^-'°-"(T5^k)(^)^ 

At this point it is instructive to compare the typical strain associated with the BWM signal from SMBH merger with the 
strain amplitude associated with th e inspiral phase prior to the SM BH merger. The characteristic strain amplitude associated 
with the inspiral phase is given by (jSathvaprakash fc Schut3l2009l ) 

for the inspiral of SMBH of comparable masses. Although the characteristic strains associated with the BWM and the inspiral 
have comparable amplitudes, they lead to different timing residual signatures. In a nutshell, the qualitative difference in the 
residuals arises because the BWM signal is characterized by a permanent offset of GW field which leads to linear growth of 
the residuals with time, whereas the inspiral signal has a quasi-periodic time varying GW field and consequently does not 
lead to linear growth of the residuals (see Section [S}. 

In order to gain further insight, it is helpful to present a heuristic argument for concentrating our effort on the signature of 
BWM signal in pulsar timing residuals instead of the inspiral signal although both have the comparable characteristic strains. 
We refer the reader to Section [3.31 for a rigorous calculation, however for the purpose of order-of-magnitude estimation, it is 
reasonable to assume that the signal-to-noise ratios (SNR) associated with PTA measurement of BWM and inspiral signal 
are related as: 

SNRbwm / -Rbwm \ 
SNR 

insp insp / 

where -Rbwm, -Rinsp are timing residuals associated with the BWM and inspiral signal, respectively. The BWM residuals 
grow linearly with time -Rbwm ~ /i™'^™Tobs (see Section [3|, where Tobs is the total duration of observation. In the case of 
an inspiral signal due to its quasi-periodic nature the residuals oscillate with a characteristic amplitude -Rinsp ~ ^'"''^/winap, 
where Winsp ~ [7.5 • 10^ (M/10**Mq) sec] is the characteristic frequency associated with the inspiral (which is the frequency 
at the Last Stable Circular Orbit). Thus, assuming Tobs = 10 yrs, the expected ratio is: 

SNRbwm h^^^Tohsi^insp „ o 1 n4 

-SNR— h^^ robsC.insp~2.10 . 

For this reason, in what follows, we shall ignore the possible contribution to the pulsar timing signal coming from inspiral of 
the SMBH binary. 



2.2 Assessment of the event rate of SMBH mergers 

The crude estimation above shows that the amplitude of BWM from SMBH mergers is within the expected reach of pulsar 
timing measurements. Thus, the prospect of detecting the BWM signal from SMBH binary mergers crucially depends on the 
expected event rate. A typical pulsar timing measurement lasts for a period of Tobs = 10 yrs. For this reason, if the event rate 
in a volume of 1 Gpc^ is larger or comparable to 0.1 events yr~^, it is likely that a pulsar timing measurement would be able 
to detect such an event. 

The event rate o f SMBH mergers remains very un certain. Various theoretical models predict rates that differ by 2-3 
orders of magnitude l| Baker et al. I I2OO7I : IVolonterillioO^ . The SMBH merger rate primarily depends on two factors both of 
which need further study. The first factor is the number of merging galaxies that contain SMBH in the mass range of our 
interest. The second factor is the fr action of gal a xy me rgers that leads to SMBH mergers. The galaxy merger rate at redshifts 
2 < 1 is around one per year (see iBaker et al.l |20o3) and references therein). Therefore, simply assuming that each of the 
merging galaxies contains an SMBH and that these SMBH coalesce in each galaxy merger, the SMBH merger rate should also 
be around one per year. Although very crude, these estimates indicate that the SMBH merger rate is in the right ball-park 
to be detected in pulsar timing measurements. 

The event rate for SMBH mergers in the mass range 10^ Mq < M < l(f Mq was calculated to be 0.4 events yr~^ in 
l|Sesana et al.ll2003 ). These mergers occur for redshifts 2 < 4, with maximum event rate at z ~ 2. In addition, their calculations 
show that at least over 20% of SMBH mergers have mass ratios larger than 0.2. A concordant number for the expected even t 
rate, of 0.1 events yr~^ for SMBH merger at redshifts z < 1 with total mass M ~ 10* Mq comes from (|Enoki et al.l 120041 '). 
Moreover, as suggested by Fig. 6b in (|Enoki et al.ll2004l '). the event rate for SMBH mergers with masses lO'^ Mq < M < 10^ Mq 
can be about 1 events yr"'^, albeit coming pr imarily from a larger redshift z ~ 3. Along with the theoretical estimates of 
the SMBH merger event rates, recent studies (|Conselice.Yang fc Bluclj|2009l ) suggest Rg ~ 10"'' events Gpc~''yr~^ for the 
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rate of major galaxy mergers. Assuming that each major galaxy merger is associated with SMBH binary coalescence, both 
numerical estimations and observational data allows us to conclude that a rate of one event per ten years of observations 
within a redshift of z < 0.5 is not unrealistic. 



3 SIGNATURE OF BURSTS WITH MEMORY IN PULSAR TIMING RESIDUAL 
3.1 Pulsar timing residual due to a BWM 

A gravitational wave propagating between the pulsar and observ e r leads to a mo dulation in the frequency of the observed 
pulsar signal given by l|Estabrook fc Wahlauistlll975l : ISazhinlll978l : lDetweiie3ll979l ) 



uo 2c J V dt 







(7) 

path 



where uq is the unperturbed pulsar frequency in the absence of gravitational waves and Ai'{t) = v{t) — vo is the variation 
of the observed pulsar frequency due to the presence of a gravitational wave characterized by the field hij. D is the distance 
from the pulsar to the observer and c is the speed of light (below we set c = 1). The expression in the brackets is evaluated 
along the light ray path from the pulsar to the observer, with the integration variable A being the distance parameter along 
this path, and e' being the spatial unit vector along the path. The unperturbed light ray path is given by 

t{\)^t-\, x\s) = xh-e'X, (8) 

where t is the time of observation and Xq is the position of the observer. Without loss of generality we set Xp = by choosing 
a spatial coordinate system with the observer at origin. 

The pulsar timing measurements primarily measure the tim ing residuals, i.e . the difference between the actual pulse arrival 
times and times predicted from a spin-down model for a pulsar (|Detweileilll979l l . The timing residual s{t), accumulated during 
a time interval of length t beginning from an initial time tin = 0, due to the presence of a gravitational wave can be calculated 
from the expression for the frequency modulation ((7|) in the following way 

s(t)=jdr^. (9) 

J 1^0 



the residual s [t) has the dimensions of time and is customarily measured in nanoseconds. 

In our analysis we shall assume that the gravitational wave source is sufficiently far from the observer on Earth in 
comparison to the distance D between the observer and the pulsar, so as to treat the GW in the plane wave approximation. 
In this case, the GW incoming from the direction given by the unit vector rii can be presented in the form 

h^j (^x\ ^ h+(t- Uix"^ p+ + ftx (t~ n^x'j pfj, (10) 

where h+ and hx are the amplitudes corresponding to two linear polarization states of a GW. The linear polarization tensors 
and pf^ can be expressed in terms of two mutually orthogonal unit vectors U and vrii lying perpendicular to the direction 
of the wave propagation Ui as follows — [hlj — rriimj) and pfj — [hmj +milj). The various projection terms that are 
encountered in evaluating expressions ((Tjl and (|9]) can be written as 

UiC^ = /I, ptj^^e^ — (1 ^ ^^)cos2(j), PijCe' = (1 — ji')sh\2(t>, (11) 

in terms of angular variables /i — cosO and (j) (see, e.g., iBaskaran erall (120081 )1. Here 9 is the angle between the direction of 
GW propagation and the direction from the pulsar to the observer e'. Angle 4> is the azimuthal angle of e' with respect to 
the principal direction U and rrii characterizing the GW polarization tensor, i.e. the angle between vector U and the projection 
of vector e* onto Zimi-plane. 

The expression for frequency modulation (O can be integrated for the GW given in the form (|10p exactly. Taking into 
account the projection factors (lll|) . we get 

= i (1 + /x) ( \h+ (t) cos 2</> + /ix (t) sin ] - 

- ^h+{t- D{1- ^^)) cos 2(p + hy, (t - D (1 - ^)) sin 20 j | . (12) 

The above expression has a clear physical interpretation, according to which the variation in frequency is directly proportional 
to the difference between the GW field strength at the place and time of observation (terms in the first square bracket) and 
its strength at the place and time of signal emission (terms in the second square brackets). In specific case of a BWM signal 
that reached the Earth during time span of pulsar observation, the second term will naturally vanish since the strength of the 
signal at the place and time of emission would have been equal to zero. On the other hand, it can be seen from ((9} and ([12} 
that if BWM front crossed the pulsar neighborhood at or before the moment of emission the second term would give rise to 
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Figure 1. The solid line shows the expected timing residual due to a BWM with /i™^'" = 5 ■ 10"^'' and = 5 yrs, in a pulsar with 
directional angles fi = <p = 0, during an observation spanning 10 yrs, after subtraction of the quadratic fitting term. The dashed and the 
dashed-dotted lines show the BWM residual signal before subtraction and quadratic subtraction terms, respectively. 



linear rising term in prefit timing residuals. This linear term would be absorbed after fitting procedure (see below) and so the 
second term can be neglected in this case also. For these reasons, below we shall omit the second term in (|12p . 

Using expression l|12|) and neglecting the terms proportional to the GW field strength at emission, we arrive at the 
expression for the timing residual 



sit) = I + f^) I i J dT h+ (t)] cos2<^ + n dr ftx (r) 



sin 20 



(13) 



We shall model the BWM in a simple analytical manner (that would be applicable in a wide range of situations) as a 
step-function signal 



h+{t) = h""'"'e{t-tB), /ix(t) = o, 

where t_g is the time the BWM signal reaches the observer on Earth. The function Q{x) is the Heaviside step function 



e{x) 



x>0 



From (|13|) and (|14p we find the prefit timing residuals: 
1 , 



^ 7 mem 

prefit 2 



(1 + ^) cos 2(j> {t - te) e (t - ta) : 



(14) 



(15) 



(16) 



Postfit residuals are obtained by the removal of linear and quadratic trends from prefit ones to take into account a priori 
unknown values of the pulsar period and its first derivative. 



/,\ ^ 7 mem / 

s{t) = -h 



'(1 + m) cos 20 I{t), 
where X (t) is given by the expression 

l{t) = (t-tB)e(t-tB) -Xquad(t). 

The quadratic subtraction term Iquad (t) can be written in a general form 

tB)+C. 



(17) 
(18) 
(19) 



Iquad (t) =a{t-tBy +b{t 

The coefficients a, b and c can be evaluated by minimizing the integral J^"^" dt 1? (t) (where Tobs is the total duration of 
observation) with respect to these variables. For example, in the case of a BWM signal occurring at ts ~ robs/2, these 
coefficients take a particularly simple form 

15 ,1 STobs 



(20) 



ler^bB 2 ' 64 

The signature in the pulsar timing residuals of such model BWM is shown in Fig. [T] 

In the following subsections, we shall assess the sensitivity of PTA observations to the BWM signal, and estimate the 
expected rate of observable events. 
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3.2 Signal-to-noise ratio 

In order to quantify the detection ability of PTA, in tiie present section we introduce and evaluate tlie signal-to-noise ratio 
for a BWM event. The observed residuals in a pulsar timing array can be presented in the form 

Rc{U) = s^{U) + nc{ti), (21) 

where the index a — 1, .., Na marks the residuals measured for the a**^ pulsar, and the index i — 1, .., Nt is the number of the 
measurement. A'^^ and Nt are the number of pulsars in the timing array and the number of time observations per individual 
pulsar, respectively. Sa (ti) is the part of residuals due to BWM (Eq. (|17|) ). and (ti) is the noise. 

We shall assume that Ua (ti) is gaussian stationary white noise, uncorrelated for different pulsars. Under these assumptions, 
the noise correlation function will have the form 



Ua (ti) np {tj) = al{a)5ijSai3, (22) 
where S^fs and Sij are the Kronecker deltas. The residuals Sa{t) can be conveniently rewritten as 

Sa{ti) = h'^^'^f (^a,(l>a)T{ti) , (23) 

where I (ti) is the common part of the signal for all pulsars (|18|) . and / (na, 4>a) is the part of the residuals that depends on 
the orientation angles of the pulsars (/!□,, </>□,) 

f {^Jia,4>a) = ^{'^+ ^J^a)cOs2(|>a. (24) 

The common way to extra ct a signal with a known form from a gaussian stationary noise is by using a matched filter 
ijSathvaprakash fc Schut3l2009l ). The signal-to-noise ratio p attainable using the matched filter can be written in terms of the 
power of the noise o"^(a) and the expected signal Sa (ti) 



0=1 V ' i=l 



p' = l.(-59-Tl.«:i(i01. (25) 



Note that the above expression for signal-to-noise ratio is derived in time-domain variables which is more convenient in our 
case. Using (|23|) . p can be rewritten in a factorized form 

(26) 



Assuming a near isotropic sky coverage for the pulsar timing array and a similar noise level for all the pulsars a„, the terms 
in the first bracket can be calculated as a sky-average value 



—y 



' ^dMd0f(M,0) = -i-. (27) 



Na ^ crn(Q) 47rcr2 J 60-2 

The terms in the second bracket in expression (|26[) can be approximated by an integral 



it). (28) 



This integral can be calculated explicitly if the form of Iquad in Eq. (|18|) is known. For example, in the case of t_g = Taha/'2, 
the form of Xquad is given by expressions (|19|l and (I20|l and 

dtl'(i) = ^, (29) 


It is convenient to introduce a dimensionless quantity t 



which is about one for < < Tobs. 

In the future PTA observati ons Nt — 250, A''c, = 20, Tobs = 10 yrs, (t„ = 100 ns (that is the level of sensitivity currently 
achieved for a number of pulsars (jVerbiest et al.ll2009l ) ). Taking these numbers along with a signal strength /i™*^™ = 10~^^ as 
guidelines and using approximations (|27|l and (I28|) . the expression for signal-to-noise ratio can be rewritten in the form 



p = 1.64 



( Nt_Y (N^y ( / 100ns \ 



10-iV\^250y V20y VlOyrsyV / 



(31) 
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3.3 Expected event rate 



Equation (|3H) shows the practical possibihty of using pulsar timing measurements to detect individual GW bursts that 
accompany SMBH mergers. As can be seen from Eq. ((31}, the BWM with an amplitude of ~ (1.5 — 2) ■ 10~^^ can be 

detected in a 10-years run of PTA observations with a signal-to-noise ratio of ~ 3. Eq. ([5]) indicates that such a signal would 
be produced by the coalescence of two SMBH with equal masses around 3.5 ■ 10* Mq at a distance of 1 Gpc. 

It is instructive to substantiate this rate by properly taking into account the contribution to event rate from SMBH 
mergers at various cosmological distances (redshifts). The rate of GW bursts of sufficient strength from coalescing SMBH 
coming from redshifts z < znm, N{ziiin), ([A^(2:iim)] = yr^^) can be calculated from the number density of SMBH mergers n{z) 
with masses Mbh > Miii^{z), {[n{z)] = Mpc~^) and the characteristic merging frequency for a SMBH, r;(z), [77] = yr"^ The 
event rate A*'(2;iini) of bursts with sufficient strength is 

NM = no(l + zf^Ai^r^dr , (32) 

Jo L + z 

where r is the metric distance determined from 

— = " ^ (33) 

dz Ho + zy-> + Qa 

in a flat ACDM cosmological model. Here .ffo is the present-day value of the Hubble parameter and Qa, stand for the 
cosmological constant and matter energy content in units of the critical density, fim + ^a ~ 1. For numerical estimations 
below we adopt the standard values Ho = 72 km s~^Mpc~^, Qm ~ 0.3 and Qa ~ 0.7. 

The local comoving dens ity of SMBH no can be estimated using the integral SMBH mass-function from the study of 
ICaramete fc BiermannI ( 20091 ) : 



(P{Mbh) = N{M > Mbh) = 7 • 10"' (^^^^^ Mpc'^ . (34) 

The local SMBH merging rate r/o can be estimated, assuming that each SMBH with the mass greater than 10* M0 has 
undergone at least one merging in the pas t, from the lo cal density of SMBH no and the specific merging rate per unit volume 
TZo- Using the common assumption (e.g. lEnoki et al.l |2004,)) that each major galaxy merging is associated with a SMBH 
coalescence, the specific SMBH merging rate can be obtained from the analysis of major mergings of galaxies with stellar 
mass M. > 10^° Mq (Conselice et al. 2009), 7^o = 10^ Gpc"^Gyr-\ 

7^o 10-^^ Mpc-^yr-^ 10 -1 .o.n 

po 10 ^ Mpc 

In this estimation we have used the local density of massive galaxies po ~ 0.01 Mpc~^. The obtained value agrees with the 
frequency of major galaxy mergers from (Conselice et al. 2009) r;o = (1 — 2) • 10"^° yr~^- 

We are interested only in the SMBH merger events that could be registered with PTA. For this reason, more distant 
SMBH mergings should be more massive to be detected, and the corresponding minimal detectable mass scales as 

Mun,{z) oc r(l + z) . 

Setting the reference point at ro = 1 Gpc [z ~ 0.2) we get the corresponding reference values for SMBH mass Mq and local 
comoving density of SMBH no: 

Afo = 3.5 • 10* Mq, no = 6 ■ 10"^ Mpc"^ (36) 
In terms of the reference values, the minimum detectable mass can now be expressed in the form: 

^-(^) = ^^"(^^[ttI!))- (3^) 

Using Eqs. ((Mil, IMJ and JST]), we arrive at 

n(,) = no ( ' = 6 . 10-^ ( ^4^^) ' - 10"' { -j^X ■ (38) 

Substituting Eq. ((3EI into Eq. (EH) we find 

iV(ziim) = 47r • 10~"rg / dr = 47r • 10""r,^ / -^dz (39) 
Jo Jo dz 



After substituting (|33[) into above equation, we finally obtain 
N(znn.) = 47r.l0-" r""^==i===d2 



Ho Jo ya^TT^pTTII 
1 



^fim{i + zy^ + nA 



dz. (40) 
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We set our l imitin g redshift to = 5, bearing in mind that SMBH assembly increasingly occurs around 2 ~ 3 (see, e.g., 
lEnoki et al] (12004 )). In that case the value of the integral in Eq. (|40p is close to 2 and the detection rate is 

iV(^iim = 5) ^ 10'=* yr-^ . (41) 

Taken at face value, this rate seems to be fairly low, but we can increase it significantly if we take into account the strong 
redshift dependence of the specific major galaxy merging rate ti{z) = rio{l + z)^ . With this factor the integral in Eq. (|4U|I 
takes the approximate value: 

1 . 



(l + z)^d2Ril.25-10"-"'*'". (42) 



This approximation slightly underestimates the integral for /3 ~ 0. The analysis of (Conselice et al. 2009) suggests /? = 2 — 3. 
Adopting /3 = 2 we obtain the SMBH detection rate from the limiting redshift = 5 

iV(2H,. = 5) ^ 6 ■ 10-4+«-58/5 yr-' ^ 10-2 yr-i . (43) 

Our figure of merit is the number of detections with SNR above a fixed threshold value p in a PTA observation with given 
sensitivity (i.e. the rms of timing residuals (t„) over the total observing time Tobs. This number is determined by the product 
N X Tobs. The detection rate N [zu^) cx no (Mo) oc Mq"^ cx (fe™'=™)"^, where /i™"™ is the signal amplitude that triggers the 
detector at a given SNR level p (see Eq. (pi) '), i.e. /i™™ oc pN^'^'^ Na^^'^T'^^^cjn. So finally, we find the number of detections 
in observations with duration Tobs is 

It can be seen that the expected number of detections is very sensitive to the duration of observations and the rms noise level. 



4 CONCLUSIONS 

We have shown that future pulsar timing measurements will be capable of detecting individual gravitational wave bursts that 
should accompany SMBH mergers. A GW burst with memory with an amplitude of ~ (1.5 — 2) • 10~^^ leaves the characteristic 
imprint in the pulsar timing residuals and can be detected in a 10-years run of Pulsar Timing Array observations with the 
present-day characteristics at a signal-to-noise ratio of ~ 3. Such a signal is expected to be produced by the coalescence of two 
SMBH with equal masses around 3.5 • 10* Mq at a distance of 1 Gpc. We estimate the rate of SMBH coalescences producing 
the BWM with such an amplitude to be around a few hundredth per year. The number of detections of such GW bursts 
from SMBH mergings by a PTA array with given characteristics is given by Eq. (|44p that shows a strong dependence on the 
total time of observations and the noise level in pulsar timing residuals. It is expected that increase in the PTA sensitivity 
in the near future will allow detection of BWM from coalescing SMBHs at a signal-to-noise ratio more than three. Future 
radiotelescopes, especially SKA, will be able to increase the sensitivity several times due to the decrease in timing residual 
noise and the extension of number of pulsars in pulsar timing program, thus making this method complementary to the LISA 
space mission for detection of coalescing SMBHs with masses > 10* Mq to which the sensitivity of LISA is reduced. 

Finally, it is worth pointing out that the analysis conducted in this paper is applicable for any source of BWM, not 
specifically restricted to GW signals from SMBH mergers. As was mentioned above, BWM are a generic feature of GW 
sources that emit energy in an asymmetric fashion. These BWM events will leave an imprint in pulsar timing residuals as 
long as they are sufficiently bright /i™^'" > 10"^^. 

When this paper was alr eady submit ted, two studies dealing with the influence of BWM on PTA appeared (|Setdl2009l : 
Ivan Haasteren fc Levinll2009l ). Results of (Ivan Haasteren fc Levinll2009l ) are based on Bayesian analysis method and give the 
same conclusions, results of \ Setoll2009h with different treatment of background noise is somewhere more pessimistic. 



ACKNOWLEDGEMENTS 

The authors thank Joseph Romano, Patrick Sutton, A.V. Zasov and O.K. Sil'chenko for useful discussions. Also the authors 
thank Kip Thorne for fruitful suggestions. The work is supported by RFBR grants 07-02-00961, 07-02-01034 and 09-02-00922. 
This research has made use of NASA's Astrophysics Data System. 



REFERENCES 

Baker J. et al., LISA Science Document LIS A-LlST-KP-436, \http://www. srl.caltech.edu/lisa/documents/lisa-Science-case.pdj 
Baskaran D., Polnarev A.G., Pshirkov M.S., Postnov K.A., 2008, Phys. Rev. D, 78, 044018 
Bertotti B., Carr B.J., Rees M.J., 1983, MNRAS, 243, 945 
Blanchet L., Damour T., 1992, Phys. Rev. D, 46, 4304 
Bontz R.J., Price R.H., 1979, Astrophys.J., 228, 560 



Observing gravitational wave bursts in pulsar timing measurements 9 



Braginsky V.B., Grishchuk L.P., 1985, Zh. Eksp. i Teor. Fiz. , 89, 744 

Braginsky V.B., Thorne K.S., 1987, Nature, 327, 123 

Caramete L.I., Biermann P.L, 2009, preprint (arXiv:09 08.2764[l 

Christodoulou D., 1991, Phys. Rev. Lett., 67, 1486 

Conselice C.J., Yang C, Bluck A.F.L., 2009, MNRAS, 374, 1956 

Cordes J.M. et al., 2004, New Astron. Rev., 48, 1413 

Cutler C, Thorne K.S., 2002, prepr int (|arXiv:gr-qc/0204090 I 
Demorest P. et al., 2009, preprint ( arXiv:0902.2968p 
Detweiler S., 1979, ApJ, 234, 1100 

Estabrook F.B., Wahlquist H.D., 1975, Gen. Rel. and Grav., 6, 439 

Enoki M., Inoue K.T., Nagashima M., Sugiyama N., 2004, Astrophys.J., 615, 19 

Favata M., 2009, Journ. of Physics: Conference Series, 154, 012043 

Foster R.S., Backer D.C., 1990, Astrophys. J., 361, 300 

Grishchuk L.P., Lipunov V.M., Postnov K.A., Prokhorov M.E., Sathyaprakash B.S., 2001, Physics-Uspekhi, 44, 1 

Hellings R.W., Downs G. S., 1983, ApJ, 265, L39 

Hobbs G., 2005, Publ. of the Astron. Soc. of AustraUa, 22, 179 

Janssen G.H. et al., 2008, in 40 YEARS OF PULSARS: Millisecond Pulsars, Magnetars and More. AIP, Conference Pro- 
ceedings, Volume 983, 633 
Jenet F.A., Lommen A., Larson S.L., Wen L., 2004, ApJ, 606, 799 
Jenet F.A., Hobbs G.B., Lee K.J., Manchester R.N., 2005, ApJ, 625, L123 
Jenet F.A. et al., 2006, ApJ, 653, 1571 

Keating B.C., Polnarev A.G., Miller N.J. and Baskaran D., 2006, Int. J. Mod. Phys. A 21, 2459. 

Kopeikin S.M., 1997, Phys. Rev. D, 56, 4455 

Kovacs S.J., Thorne K.S., 1978, Astrophys. J., 224, 62 

LIGO official website, \http://www. ligo. caltech. edu/ 

Manchester R.N., 2008, in 40^YEARS OF PULSARS: Millisecond Pulsars, Magnetars and More. AIP, Conference Proceed- 
ings, Volume 983, 584 
Nazin S.N, Postnov K.A.,1997, A & A, 317, L79 
Ott CD., 2009, Class. Quantum Gravity, 26, 063001 
Payne P.N., 1983, Phys. Rev. D, 28, 1894 

Reisswig C, Husa S., RezzoUa L., Dorband E.N., PoUney D., Seller J., 2009, preprint! arXiv:0907.0462) 

Romani R.W., 1989, "Timing of millisecond pulsar array", in Timing Neutron Stars, edited by H. Ogelman and 

E. P. J. van Heuvel, 113 
Sathyaprakash B.S., Schutz B.F., 2009, Living Rev. Relativity, 12, 2 
Sazhin M.V., 1978, SvA, 22, 36 

Sesana A., Vecchio A., Volonteri M., 2009, MNRAS, 394, 2255 

Sesana A., Haardt F., Madau P., Volonteri M., 2009, ApJ, 611, 623 

Seto N., 2009, preprint l|arXiv:0909.1379D 

Smarr L., 1977, Phys. Rev. D, 15, 2069 

Thorne K.S., 1992, Phys. Rev. D, 45, 520 

van Haasteren R., Levin Y., 2009, preprint (arXiv:0909 .0954|) 

Verbiest J.P.W et al., 2 009, preprint (arXiv:0908.0244 ) 

VIRGO official website, http://www.virgo.infn.it/ 

Volonteri M., Invited review at the Sixth International LISA Symposium, preprint I arXiv:astro-ph/060974l| 
Zel'Dovich Ya.B., Polnarev A.G., 1974, SvA, 18, 17 



